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Circuit QED on a chip has become a powerful platform for simulating complex many-body physics. In this 
report, we realize a Dicke-Ising model with an antiferromagnetic nearest-neighbor spin-spin interaction in 
circuit QED with a superconducting qubit array. We show that this system exhibits a competition between 
the collective spin-photon interaction and the antiferromagnetic nearest-neighbor spin-spin interaction, 
and then predict four quantum phases, including: a paramagnetic normal phase, an antiferromagnetic 
normal phase, a paramagnetic superradiant phase, and an antiferromagnetic superradiant phase. The 
antiferromagnetic normal phase and the antiferromagnetic superradiant phase are new phases in 
many-body quantum optics. In the antiferromagnetic superradiant phase, both the antiferromagnetic and 
superradiant orders can coexist, and thus the system possesses Z^®Z 2 symmetry. Moreover, we find an 
unconventional photon signature in this phase. In future experiments, these predicted quantum phases 
could be distinguished by detecting both the mean-photon number and the magnetization. 



Circuit quantum electrodynamics (QED) based on superconducting qubits is a fascinating topic in quantum 
optics and quantum information 1 " 3 . This artificial spin- 1/2 particle can be controlled by tuning the 
external magnetic flux and gate voltage 4 " 6 . Moreover, a strong spin-photon coupling has been achieved, 
which allows to implement quantum operations for long coherence times 7,8 . Recently, many important quantum 
effects in atomic physics and quantum optics have been observed in this artificial spin-photon interaction 910 . 
Particularly, experiments realized multiple superconducting qubits interacting with a transmission -line res- 
onator 11 . These experiments allow to explore many-body phenomena via circuit QED 12 " 21 . For example, the 
challenging Dicke quantum phase transition from a normal phase to a superradiant phase, which was predicted 
more than 30 years ago 22 " 24 , can be realized by controlling the gate voltage or external magnetic flux 25 " 28 and the 
no-go theorem arising from the Thomas-Reich-Ruhn sum rule may be overcome 29,30 .Moreover, the Jaynes- 
Cummings lattice model 31 can also be simulated by an array of transmission-line resonators, each coupled to a 
single artificial particle 32,33 . In addition, by measuring the microwave photon signature, the many-body none- 
quilibrium dynamics, as well as the known phase diagrams, could be derived 34 " 36 . 

On the other experimental side, superconducting qubits can couple with each other, forming an array with an 
effective nearest-neighbor spin- spin interaction 37,38 . Thus, it is meaningful to explore the many-body physics 
when a superconducting qubit array couples with a transmission-line resonator because there exists a competi- 
tion between the collective spin-photon interaction and the nearest-neighbor spin-spin interaction. Recently, 
sudden switchings, as well as a bistable regime between a ferromagnetic phase and a paramagnetic phase, have 
been predicted 39 , attributed to this competition. 

In this report, we investigate the quantum phases in circuit QED with a superconducting qubit array, which is 
governed by a Dicke-Ising model with an antiferromagnetic nearest-neighbor spin-spin interaction. By consider- 
ing the competition between the collective spin-photon interaction and the antiferromagnetic nearest-neighbor 
spin-spin interaction, we predict four quantum phases, including: a paramagnetic normal phase (PNP), an 
antiferromagnetic normal phase (ANP), a paramagnetic superradiant phase (PSP), and an antiferromagnetic 
superradiant phase (ASP). The ANP and the ASP are new phases in many-body quantum optics. In the ASP, both 
the antiferromagnetic and superradiant orders can coexist, and thus the system possesses Z^®Z 2 symmetry, i.e., 
both (7(1) and translation symmetries are broken simultaneously. Moreover, we find an unconventional photon 
signature in this phase which could increase from zero to a finite value and then decrease when increasing an 
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effective magnetic field. In future experiments, these predicted 
quantum phases could be identified by detecting both the mean- 
photon number and the magnetization. 

Results 

System and Hamiltonian. Figure 1 shows our proposed quantum 
network. Many superconducting qubits connected in a chain couple 
capacitively to their neighboring qubits and also interact identically with 
a one-dimensional transmission-line resonator. The corresponding 
Hamiltonian is given by 4 " 6 



Hi 



5> 



,j cos (p { , 



where w f is the number of Cooper pairs on the /th island, 



2e 



(1) 



(2) 



is the dimensionless gate charge with gate capacitance C g and gate 
voltage V ff 



Ej = 2E°j cos ( n 



(3) 



is the tunable Josephson tunneling energy with single-junction 
Josephson energy Ep external magnetic flux and magnetic flux 
quantum <J> 0 , (pi is the phase of the superconducting order parameter 
of the /th island, and C is the capacitance matrix. The element C u = C E 
is the total capacitance connected to the ith island and C u ± i = — C is 
the coupling capacitance between two adjacent superconducting qubits. 
In general, the coupling capacitance C is much smaller than the total 
capacitance C E . As a consequence, the next-nearest-neighbor term of C 
can be neglected safely and the Hamiltonian Hi is rewritten as 



2e 2 
— ( 



~ N g) ~Ej cos (Pi 



4e 2 C 

C 2 



5>,-N g )(n, + I -^.(4) 



Near the degeneracy point, only a pair of adjacent charge states («,■ = 0 
and rii= 1) on the island are relevant. If we define these charge states as 
the effective spin basis states | j,) and i.e., \n { = 0) <-> | | z ) and |«,- = 
1) \li)> me Hamiltonian H 2 reduces to the form 



H 3 

where |T,XT/I 
spin operators, 



=/£^ +1 +*£-i-y£<4 



■\U)(U\ and 0 i=|T l >a i |-U,><t,| arethePauli 



N„ 



is an effective magnetic field, and 



/ = 



(6) 



(7) 



describes the capacitance-induced nearest-neighbor spin-spin 
interaction. 

Now a one-dimensional transmission-line resonator is placed in 
parallel to the superconducting qubits. All superconducting qubits 
are situated at the antinode of the magnetic field induced by the 
oscillating supercurrent in the transmission-line resonator 2 ' 3 . Due 
to the boundary condition at the end of the transmission-line res- 
onator, these superconducting qubits are controlled only by the mag- 
netic component, which shifts the original magnetic flux O x by 



where 



71 V ' 



TlSo ihlCD 

d<D 0 V L 0 ' 



(8) 



(9) 



/ is the inductance per unit length, and S 0 is the enclosed area of the 



superconducting qubit. In the Lamb-Dicke limit (go y {erf a) + l^Cl), 
together with the condition O x = <D 0 /2, an effective Hamiltonian for 
Fig. 1 is obtained by 



tf=/£<< +1 + S £<+^£<(a + «t)W, 

i i v ^ i 



where 



(10) 



(11) 



is the collective spin-photon coupling strength. In this Hamiltonian, 
all parameters can be controlled independently. For example, the 
effective magnetic field s can be tuned via the gate voltage V g from 
the negative to the positive. For simplicity, we address mainly the 
case 8 > 0 in the following discussions. 

The Hamiltonian (10) is a Dicke-Ising model with an antiferro- 
magnetic nearest-neighbor spin-spin interaction. This Hamiltonian 
shows clearly that when both g and / coexist, the collective spin- 
photon interaction has a competition with the nearest-neighbor 
spin-spin interaction. As a result, it exhibits exotic phase transitions 
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beyond the previous predictions of the standard Dicke (Ising) model. 
For example, a first-order superradiant phase transition has been 
predicted 40 ' 41 , when / < 0. In this report, we will find rich quantum 
phases including the PNP, the ANP, the PSP, and the ASP for / > 0. 

Quantum phases. For the Hamiltonian (10), the quantum phases 
can be revealed by calculating the ground-state energy and the order 
parameters via a mean-field approach 42 . In the classical picture, the 
spin in the Hamiltonian (10) can be represented as a vector line in the 
xz plane with the unit vector n = (cos (pfi, sin (p). Thus, we can 
introduce a variational ground-state wave function 



where 



and 



|W>=n|^>®|^>, 



w=MM).-(j-!)] ! 



(12) 



(13) 



|A 0 >=exp(-|+A 0 /)|0) (14) 

are the spin and boson coherent states, to describe both the 
antiferromagnetic and superradiant properties. 

Since the antiferromagnetic exchange interaction (/ > 0) leads to a 
staggered arrangement of all spins in the z direction, we should 
consider two sublattices with q> 1 and — (p 2 , which corresponds to 
the odd and even sites of spins, respectively, in the ground-state wave 
function. After a straightforward calculation, the scaled ground- state 
energy 



E={k 0 ,(p 1 ,(p 2 \H\k 0 ,(p 1 ,(p 2 )/N (15) 



is given by 



E = wk 2 — /sin (p x sin + 

e ( 16 ) 
- (sin ^ — sin<^ 2 ) +g2(cos<^ 1 + cos#? 2 ), 

where k = k Q j Vn and the parameters (k, <p x and q> 2 ) are to be 
determined. 

As shown in the Methods section, by minimizing the ground-state 
energy E with respect to the variational parameters (k, <p iy <p 2 ), we 
obtain three equilibrium equations: 

wk+g(cos a cos/?) =0, (17) 



/(sin 2a - sin 2fS) - s cos(a + ft) - 2gk sin(a + ft) = 0, (18) 



/(sin 2a - sin 2/?) + e cos(a -fi)- 2gk sin(a -fi) = 0, (19) 



where 



71 7V 

2 '2 



71 71' 

2'2 



(20) 



(21) 



These equilibrium equations, together with the stable conditions (see 
the Methods section), determine the ground-state energy in Eq. (16) 
and the order parameters, such as the mean-photon number (aft a), 
the magnetization (S z ) and the staggered magnetization (M 5 ), which 
are given respectively by 



( flTfl ) £ 2 2 2 o 

= f^cos z acos z ^ 

%^ = ME/<)= cosasin ^ 

(M 4 ) = i|E < (-l) < (af)| = |sinacosj5| 



(22) 



The introduction of the order parameter, the staggered magnetiza- 
tion (M 5 ), is to conveniently discuss the antiferromagnetic properties 
of the Hamiltonian (10). After the ground-state energy, and espe- 
cially, the order parameters, are obtained, several rich phase dia- 
grams can be obtained. 

We first address two known limits. The first is the case when / = 0, 
in which the Hamiltonian (10) reduces to the standard Dicke model 43 



H D = £^4+-^^4(a + a t )+coa t a. 



(23) 



By means of the equilibrium equations (17)-(19) and the stable 
conditions, we find 



oc = 0,fi = 

for g<g c and 

a = 0,/? = — arcsin 
for g > g c , where g c = ^£(o/2, i.e., 



2(J+g 2 /co) 



(eft a) 
N 



= 0. 



N 



= -1 



for g<g c and 



(aft a) 
N 



4g 2 



' N 



2g 2 



(24) 



(25) 



(26) 



(27) 



for g > g c . This means that a second-order quantum phase transition 
from the normal phase (g < g c ) to the superradiant phase (g > g c ) 
occurs 44 ' 45 , as shown in Fig. 2(a). Moreover, the Dicke model has (7(1) 
symmetry in the normal phase. Whereas, in the superradiant phase 
the system acquires macroscopic collective excitations governed 
mainly by the collective spin-photon interaction term g^y^jG^ 
(a-\-cft), and thus it has Z Z 2 ®Z 2 symmetry, where Z z 2 is the global 
rotation of n around the z axis 46 and z 2 is the change of sign of the 
boson coherent state (\k 0 ) — > —\k 0 )). In experiments, this quantum 
phase transition has been observed 47 " 49 in an optical cavity with a 
Bose-Einstein condensate by measuring the mean-photon number 
(aft a) and the magnetization (S z )/N. Recently, it has been well inves- 
tigated in many-body circuit QED 25 " 30 and spin- orbit- driven Bose- 
Einstein condensate 50 . 

Forg = 0, the Hamiltonian (10) turns into the Ising model 51 



(28) 



in which a first- order phase transition from the paramagnetic phase 
to the antiferromagnetic phase at the critical point J c = e/2 can be 
recovered. In the paramagnetic phase (/ < / c ), the ground-state wave 
function is |... |||||||| ...), which implies that the system has 
translation symmetry and (S z )/N = — 1 and (M s ) = 0. In the anti- 
ferromagnetic phase (/ > / c ), the ground- state wave function 
becomes | . . . TitlTiTI • • ■)> m which translation symmetry is broken 
and (S z )/N= 0 and (M s ) = 1. 

If both g and /are non-zeros, we find four different regions: (i) a = 
0, P = -tt/2, (ii) a = +7i/2, 0 = 0, (iii) a = 0, ft ^ 0, and (iv) a ^ 0, ft 

0. Specially, the order parameters in these four regions are given 
respectively by 
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(in)- 



(a ] a) 



N co w 



, s W a ) (S z ) , , 



N N x 1 



' N 



(29) 



(30) 



,{M S ) = 0, (31) 



(v) 



(a* a) 



^ - cos 2 a cos 2 /i, = cos a sin /?, {M s )=\ sin a cos /? | , (32) 



where 



- e sin /? 



2(J-g 2 cos 2 P/w)' J 



2 o W 

cos z P= — 



7-|wi- 



CO/ 



(33) 



(34) 



In terms of the different properties of the order parameters, the cases 
(i)-(iv) are denoted by PNP, ANP, PSP, and ASP, respectively. The 
ANP ((a f a)IN = 0, (S z )/N = 0 and (M s ) = 1) and the ASP ((a f a)/N ¥= 0, 
(S z )/N 7^—1 and (M s ) 1) are new phases in many-body quantum 
optics. In Fig. 2, we plot phase diagrams for different antiferromagnetic 
spin-spin nearest-neighbor interactions. This figure shows clearly that 
these predicted quantum phases can be driven by the collective spin- 
photon coupling strength g, the antiferromagnetic nearest-neighbor 
spin-spin interaction /, and the effective magnetic field s. Especially, 
the region of the ASP becomes larger when increasing /. 



In Fig. 3, we plot phase diagrams as functions of the antiferromag- 
netic nearest-neighbor spin-spin interaction /and the collective spin- 
photon coupling strength g for different effective magnetic fields (a) 
8 = 0 and (b) 8 = col A. In the absence of s, Eq. (10) reduces to the form 



+ ) +coa ] a 



(35) 



In such a case, only the ANP and the PSP can be found, as shown in 
Fig. 3(a). When increasing g, four quantum phases are predicted 
again, as shown in Fig. 3(b). In addition, by means of the ground- 
state energy, we find that all transitions between these different 
quantum phases in Figs. 2 and 3 are of second order. 

Symmetry. In order to better understand these predicted quantum 
phases, it is necessary to discuss the corresponding symmetries. For 
the PNP and the PSP, the system properties are similar to those of the 
normal phase and the superradiant phase in the standard Dicke 
model, i.e., the system displays both (7(1) and translation 
symmetries in the PNP, and becomes Z|®Z 2 and translation 
symmetries in the PSP. However, in the ANP, though no photon is 
excited, the antiferromagnetic order emerges. This implies that in 
such a case only (7(1) symmetry can be found. Interestingly, in the 
ASP the Hamiltonian (10) is governed mainly by the term 



Hp, 



(36) 



in which there exists a competition between the antiferromagnetic 
nearest-neighbor spin-spin interaction and the collective spin- 
photon interaction. As a result, both the antiferromagnetic and 
superradiant orders coexist and the system possesses Z^0Z 2 
symmetry, i.e., both (7(1) and translation symmetries are broken 
simultaneously. 



Possible experimental observation. We first estimate the 
parameters for experiments. When we choose C £ ~ 600 aF, C ~ 
20 aF, E°j ~2n x 3.5 GHz, S 0 ~ 1 fim 2 , d ~ 10 ^m, L 0 ~ 19 mm, 
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Figure 3 | Phase diagrams for different effective magnetic fields. The plotted parameters are chosen as (a) s/co = 0 and (b) s/co = 0.25, respectively. In 
(a), the phase boundary is determined by J c = f?/co. In (b), the phase boundaries (from top to bottom) are given respectively by 

Jo'- ' ' ' " 



■ (g 2 /u+Vz 2 +g 4 /™ 2 ) /2> {Jc +g*M -g 2 /\coJ c ) = e/2, and/ c = 8/2 - g/co. 



ft) ~ 27i X 6.729 GHz, and N = 100 11 , the antiferromagnetic nearest- 
neighbor spin-spin interaction parameter and the collective spin- 
photon coupling strength are given respectively by / ~ In X 
2.2 GHz and g ~ In X 1.5 GHz (g 0 = 0.01 is responsible for the 
Lamb-Dicke approximation). In addition, the effective magnetic 
field s can range from 0 to In X 6.8 GHz by controlling the gate 
voltage V g . These parameters ensure that the system should probe the 
predicted phase transitions. To observe these phase transitions, the 
relaxation time 7\ and the coherence time T 2 should be much smaller 
than the lifetime 1/k of the photon, i.e., 7\ > T 2 > II k = 23 A ns, 
where k is the decay rate of the photon. This restriction can be easy to 
satisfy in current experimental setups (for example, 7\ = 7.3 /us and 
T 2 = 500 ns in Ref. 52). 

We now illustrate how to identify these different quantum phases. 
Here we propose to detect four phases by measuring both the mean- 
photon number (a* a) and the magnetization (S z ). For the PNP and 
the ANP, we can separate these by directly observing the magnetiza- 
tion because (S z )/N = - 1 in the PNP and (S z )/N = 0 in the ANP, as 
shown in Fig. 4(a). For the PSP and the ASP, both the photon and the 
spin are collectively excited. Moreover, when increasing the effective 
magnetic field e, (S z )/N always decreases. This means that it is dif- 
ficult to distinguish the PSP and the ASP by measuring (S z )/N. 
Fortunately, we find that in the ASP the mean-photon number has 
an unconventional behavior that could increase it from zero to a 
finite value and then decrease, as shown in Fig. 4(b). The relevant 
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Figure 4 | Order parameters versus the effective magnetic field. In (a) and 

(b), the magnetization {S^/N and the mean-photon number (a!d)/Na& 
functions of the effective magnetic field e are plotted, respectively, when 
Jlco = 0.30 and glco = 0.25. 



physics can be understood as follows. When the effective magnetic 
field 8 is applied, it can initially promote the arrangement of all spins 
from the antiferromagnetic to the paramagnetic terms 53 . For 
example, in the case of a weak e, the spin arrangement becomes 
TIT UIH TIT •••) fro m the antiferromagnetic case |... 
TITITITI •••)• This process is helpful for achieving photon-induced 
collective excitations. Thus, the mean-photon number can increase. 
However, the rearrangement of spins gives rise to an opposite result 
of the magnetization, i.e., it decreases when increasing e. For strong £, 
this effective magnetic field in the z axis leads to a large spin imbal- 
ance | . . . TTIIIIII •••) an d thus suppresses the spin-photon collect- 
ive excitations, i.e., both {a 1 a) I 'N and (S z )/N decrease when increasing 
8. In terms of the different behaviors of both (a* a) IN and (S z )/N, we 
argue that our predicted quantum phases can be identified. 

It should be pointed out that the microwave photon in supercon- 
ducting circuits is not easy to measure directly 54 via photon-number 
detectors, because its energy (Hat) is very small. However, in the 
dispersive region A ^> g, where A = 8 — co, the photon number 
can be detected by the photon-number-dependent light shift (the 
Stark plus Lamb shifts) of the atom transition frequency 55 . 
Unfortunately, to achieve the predicted phase diagrams, the system 
should work at the quasi-dispersive-strong region A/g < 4. In such a 
region, the above approach to detect photons does not work. 
Recently, it has been proposed 56 ' 57 to detect the photon by irreversible 
absorption of photons. In these proposals 56,57 , the absorbers along the 
waveguide are built with bistable quantum circuits, and can produce 
a large voltage pulse when the photon decays into a stable state. This 
suggests that in future experiments the mean-photon number could 
be detected, and then our predicted phase diagrams could also be 
observed. 



Discussion 

Let us here address the no-go theorem in quantum optics. This 
no-go theorem, demonstrated 58 in 1975, shows that in a typical 
optical cavity with an ensemble of natural two -level atoms, the 
phase transition from the normal phase to the superradiant phase 
is forbidden by the A 2 term, where A is the vector potential. 
Recently, the no-go theorem has been addressed 29 ' 30 in circuit 
QED, with many superconducting qubits interacting with a quan- 
tized voltage (microwave photon). However, in this report, the 
required microwave photon is generated from the quantization 
of the magnetic flux. In such a case, no A 2 term can occur, i.e., 
the no-go theorem is not valid. 

We now consider how the decay of the photon and the disorder in 
fabrication affect the predicted phase diagrams. When considering 
the decay of the photon, the stationary mean-photon number, which 
can be derived from 
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the presence of a strong external magnetic field (i.e., for a weak 
magnetic field, the disordered phase can occur). The phase dia- 
grams predicted here should be observable under a strong mag- 
netic field (See Figs. 2(c) and 2(d)). This means that these 
predicted quantum phases will not be qualitatively affected by 
weak disorder. 

Mean-field predictions become more accurate for larger number 
of spins. However, mean-field is often a good starting point, and 
provides some basic insight in the system. Moreover, for current 
experimental techniques, the spin number is not sufficiently large, 
but this should change in the future. For smaller number of spins, we 
can perform direct numerical diagonalization to discuss the ground- 
state properties. In Fig. 5, we plot the order parameters (a f a)/N, (S z )/ 
N, and (M s ) as functions of the collective spin-photon coupling 
strength g and the effective magnetic field s. This result shows clearly 
that for a small number of spins, the predicted quantum phases still 
exist, but the phase boundaries are affected significantly. 

In summary, we have investigated the Dicke-Ising model with 
an antiferromagnetic nearest-neighbor spin- spin interaction in 
circuit QED for a superconducting qubit array and predicted four 
quantum phases, including the PNP, the ANP, the PSP and the 
ASP, with different symmetries. Moreover, all transitions between 
these different quantum phases are of second order. We have also 
found an unconventional photon signature in the ASP, where 
both the antiferromagnetic and superradiant orders coexist. We 
believe that this system allows to explore exotic many-body phys- 
ics in quantum optics and condensed-matter physics because it 
has an interesting competition between the collective spin-photon 
interaction and the nearest-neighbor spin-spin interaction. 



Figure 5 | Order parameters via the direct numerical diagonalization 
method. In (a)-(c), the mean-photon number (d ] 'a)/N y the magnetization 
{S z )/N, and the the staggered magnetization (M 5 ) as functions of the 
collective spin-photon coupling strength g and the effective magnetic field 
s are plotted, respectively, when J/co = 0.1 and N = 7. 



da 



= coa H — — S x - iKa - 



becomes 



(37) 



Methods 

Three equilibrium equations and the corresponding stable conditions. Here we 
present detailed calculations on how to derive the three equilibrium equations (17)- 
(19) and the corresponding stable conditions. After minimizing the ground- state 
energy E in Eq. (16) with respect to the variational parameters (/I, (pi and (p 2 ), we 
obtain three equations: 



— =20JA+g{COS(p l + cos(p 2 ) = 



o, 



cE 



-] cos(p Y sm(p 2 + -cos^ — glsm(p l =0, 



(41) 



(42) 



N(oj 2 + K 2 



:(s x ) 2 



(38) 



whereas the stationary value of (S x ) remains unchanged (i.e., it is 
identical to the case without photon decay). This means that we 
can use an effective Hamiltonian 

Heff =/ V crly- 1 + -f = V o* x [a + a*) + co eff a V (39) 



with 



COeff = Vo) 2 + K 2 , 



(40) 



to discuss the phase diagrams 49 induced by the decay of the photon. 
Since the decay rate k of the photon (~MHz) is far smaller than the 
other parameters (~GHz), the decay of the photon has almost no 
effect on the predicted phase diagrams. 

In addition, the imperfections in fabrication result in a weak 
randomness in the antiferromagnetic nearest-neighbor spin-spin 
interaction 14 . Moreover, the disorder antiferromagnetic interaction 
generates a disordered phase, such as a random singlet phase 59 . In 
this phase, most spins form a singlet pair with nearby spins, and 
the residual induce weak long-distance pairs. Unfortunately, the 
disordered phase is only a local correlation, and is thus unstable in 



dE 8 
- — = — / sin (p x cos (p 2 + - cos (p 2 —gA sin (p 2 = 0. 

C(p 2 2 



(43) 



After defining a = {(p x + q> 2 )/2 e [ — tt/2, 7c/2] and/3 = {(p x - q> 2 )l2 e [-7c/2,7c/2],the 
three equilibrium equations (17)— (19) are easily derived. 

In addition, by means of three parameters {X, a and ft), the ground- state stability is 
determined by the 3X3 Hessian matrix 



M~- 



" d 2 E 


d 2 E 


d 2 E~ 




dldoi 




d 2 E 


d 2 E 


d 2 E 


deed! 


dot 2 




d 2 E 


d 2 E 


d 2 E 


_dfidk 2 


dfJda 2 





(44) 



If M is positive definite (i.e., all eigenvalues h { of M are positive), the system is located 
at the stable phases. If M is negative definite (i.e., all eigenvalues hi of M are 
negative), the system is unstable. 
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